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ABSTRACT 

The spin temperature of neutral hydrogen, which determines the optical depth and bright- 
ness of the 21 cm line, is determined by the competition between radiative and collisional 
processes. Here we examine the role of proton-hydrogen collisions in setting the spin tem- 
perature. We use recent fully quantum-mechanical calculations of the relevant cross sections, 
which allow us to present accurate results over the entire physically relevant temperature range 
1-10'' K. For kinetic temperatures Tk ^ 100 K, the proton-hydrogen rate coefficient exceeds 
that for hydrogen-hydrogen collisions by about a factor of two. However, at low temperatures 
(Tk ^ 5 K) H-H+ collisions become several thousand times more efficient than H-H and 
even more important than H-e^ collisions. 
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1 INTRODUCTION 

The 21 cm transition is potentially a powerful probe of the pre- 
reionization intergalactic medium (IGM) because of the eno rmous 
amount of neutral hydrogen in the Univer se at that time jField 
ll958l : IScott & Reeslll99(ll : ]Madau et aljl997h . It can teach us about 
reionization, the formation of the first structures and the first 
;alaxies, and even the "dark ages" before these objects formed 
Furlanetto et alJlzOOd and references therein). It is therefore cru- 
cial to understand the fundamental physics underlying the 21 cm 
transition. One critical aspect is the spin temperature, which is de- 
termined by the competition between the scattering of cosmic mi- 
crow ave background (CMB) photon s, the scattering of Lya pho- 
tons jWouthuvsenll 19521 : iFieldl 1958h . and collisions. When CMB 
scattering dominates, the IGM remains invisible because the spin 
temperature approaches that of the CMB (which is used as a back- 
light). 

Before star formation commences, collisions are the only 
way to break this degeneracy. The total coupling rate is deter- 
mined by collisions with other hydrogen atoms, protons, and elec- 
trons. At the low residua l electron fractio n expected after cosmo- 
logical recombination tSeager et alj|l999h . H-H collisions domi- 
nate. Spin exchange in such i nteractions has received a great deal 
of attention over the years (Purcell & Field' 'l956"; 'Smith' 'l96^; 
Allison & Dalgarno 1969; Zyaelman 2005; Hirata & Sigurdson 
2006) . We have also recently re-examined spin-exchange in H- 



e~ collisions using acc urate numerical cross-sections including the 
L — 0-3 partial waves dFurlanetto & Furlanettol 20071) and showed 
that such collisions become important when the ionized fraction 
a;i>0.01. 
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However, proton-hyd rogen collisi ons have received little at- 
tention in the literature; iSmithI j 19661 ) provided the most recent 
evaluation of their spin-exchange rate coefficient, but he used only 
semi-classical estimates of the cross sections. In the intervening 
years, atomic physicists have calculated the relevant quantum me- 
chanical cross sections to high accura cy using increasingly so - 
phisticated numerical techniques (e.g .. Hunter & Kurivan 197 
Hodges & Breig 199li; lKrstic & Schul tz 1999b; Krstic et al. 20O 
Our purpose here wil l be to genera te similarly accurate rate coeffi- 
cients from the Krstic et alj j2004) cross sections for use in 21 cm 
calculations. We will show that, at sufficiently low temperatures, 
H-H^ collisions could dominate the spin coupling, but that in more 
realistic circumstances they provide only a small correction to the 
usual calculation. 

The remainder of this paper is organized as follows. In Sj2l we 
briefly review the 21 cm transition. Our main results are contained 
in SjS] where we calculate the spin de-excitation rates for H-H+ 
collisions. We conclude and discuss applications to the 21 cm tran- 
sition in the high-redshift IGM in 21 



2 THE 21 CM TRANSITION 

We review the relevant characte ristics of the 21 cm tran sition here; 
we refer the interested reader to lFurlanetto et alj j2006h for a more 
comprehensive discussion. The 21 cm brightness temperature (rel- 
ative to the CMB) of a patch of the IGM is 
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where 5 is the fractional overdensity, xm is the neutral fraction, Ts 
is the spin temperature, is the CMB temperature, and dii||/dr|| 
is the gradient of the proper velocity along the hne of sight. The 
last factor accounts for redshi ft-space distortions jBharadwai & Alii 
l2004lBarkana & Loebll2005h . 

The spin temperature Ts is determined by competition be- 
tween scattering o f CM B photons, scattering of UV photons 
dWouthuvsenl 1 19521 : iField 1958) , and coUisions l iPurcell & Fieldi 
Il956t) . In equilibrium. 



^ ^ 1+Xc+Xa 



(2) 



where Tk is the kinetic temperature of the gas, Tc is the color tem- 
perature of the radiation field at the Lya transition, and the Xi are 
coupling coefficients. The last part of this equation describes the 
Wouthuysen-Field effect, in which the ab sorption and re-emission 
of Ly g photons mixes the hype rfine states I Wouthuysen 1952|; lFieldl 
1 19581) . We refer the reader to iFurlanetto et al. (2006.) , and refer- 
ences therein, for more detail on this component. 

The factor Xc is the total collisional coupling coefficient, 
including H-H, H-e~, and H-H'*' collisions. In this paper, we 
will focus only on the contribution from proton-hydrogen col- 
lisions , which we will denote a::^^; the o ther components are 
X™ and x'^^, with obvious meanings (see IZvgelmanll200^ and 
IFurlanetto & Furlanetto! [20071 for the most recent estimates). In 
more detail. 



(3) 



where is the spin de-excitation rate in proton-hydrogen colli- 
sions, Up is the local proton density, T, = hi^2i/kB = 0.068 K, 
i'2i is the frequency of the 21 cm line, and Aio = 2.85 x 10~^^ s~^ 
is the Einstein-A coefficient for that transition. 



3 PROTON-HYDROGEN COLLISIONS 
3.1 The spin exchange cross section 

We must first compute the cross section for spin exchange as a func- 
tion of collision energy. To do so, we begin by noting that the radial 
wave functions in the H-H+ system satisfy the Schrodinger equa- 
tions in uncoupled partial waves of angular momentum L. In the 
energy range with which we are concerned, these become 



L(L + 1) 
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where R is the internuclear distance, ^ is the reduced mass, Ei 
is the adiabatic electronic potential for the Isag (gerade) or 2pau 
(ungerade) states of HJ, and ^ ^ is the radial wave function for the 
appropriate channel. 

We assume the process to be electronically elastic, so that the 
scattering problem can be reduced i n the usual way to the compu- 
tation of phase shifts 5£ and 5^ (see IFurlanetto & Furlanettoll2007l 
for a description of the analogous transformation for H-e~ scatter- 
ing). The so -called charge tran sfer cross section act may then be 
written (e.g.. lKrstic et al.ll2004) 



(5) 



where k is the center-of-mass momentum. This cross section is de- 
fined with reference to an experiment in which a beam of polarized 
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Figure 1. Charge transfer cross section (in units of ttOq) for H-H+ coUi- 
sions, as a function of the collision energy. The solid curve shows Cct from 
iKistic et al] i2004 . The dashed curve in the inset shows the effective ranee 
approximation for the low energy behavior (see text). 



protons (with, say, spin -1-1/2) is incident on an unpolarized col- 
lection of hydrogen atoms and in which the spins of the scattered 
protons are measured, so that protons with spin -1/2 can be unam- 
biguously determined to have originated in the target. It tends to the 
usual charge transfer cross section in the classically distinguishable 
particle Umit jKrstic & Schultzill999aD . 

This cross section must be calculat ed numerically, and there 
is a long history of such attempts faunter & KurivanI 1 19771 : 
iHodges & B"rei3ll99ll : lKrstic & Schultzj |l999bV For 10"^ eV < 
E < 10^ eV, we use the most recent calculations, from lKrstic et al.l 
|2004 which are accurate numerically to six significant figures, al- 
though the neglect of coupling to higher electronic states reduces 
the accuracy somewhat at the upper end of the energy range. 

Figure [T] shows the resulting cross section, in units of Trap 
(where ao is the Bohr radius) over this entire energy range. Note 
the rich structure in act, especially at moderately small energies. 
This is in sharp contrast to the H-e~ cross section, which is smooth 
until the neighborhood of the n = 2 excitation threshold (see Fig. 
1 in IFurlanetto & Furlanettdl2007h . Capturing this detailed struc- 
ture required the use of a closely-spaced energy grid (664 points 
were used here) and the i nclusion of high-L partial waves (up to 
Lniax = 3200 at 100 eV: lKrstic et"ai]|2CI04t) . The narrow features 
are caused by shape resonances in the electronic potentials and have 
the usual Fano line shape. The broader features, known as "Regge 
oscillations," have a more subtle origin. They are generated by com- 
binations of one to three poles in the S'-matrix, each of which corre- 
sponds to one of the L = vibrational bound states of H J . These 
features are difficult to understand in the traditional partial wave 
representation, but their nature has recently been e xplained using 
the Mulholland representation jMacek et"al]|2004l) . Note that the 
charge transfer cross section shows significantly less structure than 
the elastic cross section, which also includes broad oscillations at 
higher energies due to the glory effect (Child 1984). 
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Figur e 2. Spin-exchange rate coefficients, (a): Rate coeffic i ents fo r H-H+ collisions (solid curve), H-e+ collisions (dotted curve, from lFurlanetto & Furlanettol 
120071) . and H-H collisions (dashed curve, from lZveelmM] )2005h ). (by. Ratios between the three rate coefficients in panel a. 



These calculations ignored non-adiabatic and relativistic ef- 
fects, so they were not extended below ~ 10"* eV, where such 
phenomena become significant. Instead, we use the numericall y 
converged low-energy cross sections of iGlassgold et alj 1 200f ). 
These calculations used the same technique as Krstic et all i 2004 ), 
and thus also neglected non-adiabatic and relativistic effects, but 
they at least provide us with a baseline estimate of the collision 
behavior at small energies. They show that only the L = terms 
remain non-zero at low energies, and 5q nyr as £ ^ 0. In this 
limit. 
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which is one-quarter of the total elastic cross section. Meanwhile, 
effective-range theory yields an excellent approximation for 5q: 



A;cot(5o = \ ^-fcH k In 

a 3a^ 3a 



(7) 



where the scattering length a — 801.2 ao and the polarizability 
Q = 4.5 a n paramet erize the polarization potential of the hydrogen 
atom Joiassgold et alj2005l) . 

The resulting cross section is shown for 10"^" eV < E < 
10~^ eV by the inset in Figure [T] We use the effective-range ap- 
proximation for E < 2x 10^^ eV and linearly inte rpolate the cross 
sectio n between that value and 10^* eV (where IGlassgold et al.l 
l2005l showed that the effective range approximation for the total 
cross section breaks down). This effective range theory approxima- 
tion provides a good match to numerical calculations of the spin- 
exchange cross section (using methods similar to Glassgold et al. 
l2005h in this energy regime (P. Krstic, private communication). 
Although we do not know how strongly the relativistic and non- 
adiabatic corrections will affect the cross sections, we do not expect 
them to have a dramatic effects on our final results, as discussed be- 
low. 



3.2 The spin-exchange rate coefficient 

To apply this to H-H^ collisions in the IGM, we must thermally av- 
erage the above cross sections. The spin-exchange rate coefficient 
relevant for our purposes is 



^pH 



(8) 



where fcs is Boltzmann's constant, Tk is the kinetic temperature 
of the IGM, and the prefactor with the square root is the thermally- 
averaged velocity. The factor 3/16 appears because the cross sec- 
tion as defined above assumes a polarized beam of electrons, while 
here we care about the net rate at which spin de-excitation occurs 
(Smith 1966; Krstic & Schultz 1999a). We must therefore multiply 
(Tct by 1/4 (the probability that the target atom has F — 0) and 
then 3/4 (the probability that the scattered atom has F = 1). 

The solid line in Figure |2}; shows over the full temper- 
ature range of interest. It is much smoother than act, because the 
thermal averaging washes out the resonances. At moderate to high 

1/2 

temperatures, the rate coefficient is roughly proportional to Tj^ , 
because act is itself falling only slowly with energy. However, at 
lower temperatures the cross section increases rapidly, so that 
is actually the largest rate coefficient at Tk J; 3 K. This is in sharp 
contrast to hydrogen-hydrogen and hydrogen-electron collisions, 
whose rate coefficients are shown by the dashed and dotted lines, 
respectively. The H-e~ cross section increases as E — > 0, but only 
modestly (see Fig. 1 of IFurlanetto & Furlanetto |2007|) , while the 
H-H rate decreases rapidly as iJ ^ (see below). 

Figure [2}) compares these three processes in more detail, 
showing the ratios between the various cross sections. At extremely 
small temperatures, H-H^ collisions can be the dominant pro- 
cess. Electron-hydrogen collisions quickly come to dominate at 
Tk ^ 3 K, and remain more important at all higher temperatures. 

Naively, we would expect that Ki^/k'^^ ~ y^mn /^nie ~ 30. 
In reality, the ratio is several times smaller than that, because 
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Table 1. Proton-hydrogen spin de-excitation rate coefficients 
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proton-hydrogen collisions have a substantially larger cross section 
throughout this temperature rangeQ 

Figure [21) also shows that, at high temperatures, H-H"*' colli- 
sions are marginally more efficient than H-H collisions (by about a 
factor of two); the shapes of these two rate coefficients mirror each 
other quite closely in this regime. However, at lower temperatures, 
proton-hydrogen collisions are vastly more efficient. As we have 
seen, the L = term in a^t approaches a large constant value; 
in the analogous H-H problem, the corresponding quantity actu- 
ally approaches zero be cause of an accidental cancellation in the 
5- wa ve cross sections dZvaelmM] |2005| ; ISigurdson & Furlanettd 
1200^) 

For use in other calculations, we present in Table[T]our results 
for as a function of temperature. These are fully numerically 
conve rged to the six significant figure accuracy of the lKrstic et al.l 
( 12004 ) cross sections, but as emphasized above they do not include 
all the relevant physical processes across the entire temperature 
range. In particular, at the lower end of the temperature range, non- 
adiabatic processes and relativistic corrections lead to uncertainties 
at the per cent level. For example, assuming that the cross section 
is constant below lO"* eV changes our final results by (5, 1, 0.2) 
per cent at Tk ~ (1, 3.7, 10) K, respectively. At the upper end of 
the temperature range, rotational coupling of the 2pau and 2pnu 
states affects the rate coefficients at the ~ 0.1 per cent level. 



4 DISCUSSION 

We have computed the rate coefficients for spin-exchange in H-H+ 
collisions using recent, fully quantum mechanical solutions for the 
relevant cross sec tions over most of the energy range of interest 
jKrstic et alj2004l) . We have also extended the calculation to lower 
te mperatures (Tk 1 K ) by using the approximate cross sections 
of lGlassgoldetalJj2005l) ; while these ignore non-adiabatic effects 
and relativistic corrections, they appear to be accurate to several per 
cent. 

Our results (collected in numerical form in Table [TJ will be 
particularly useful for calculating the spin temperature of the high- 
redshift IGM before reionization (and hence brightness in the 21 
cm line). Before the first luminous sources appear, this is deter- 
mined purely by the competition between collisions and CMB scat- 



tering. In most cases, protons contribute a small but non-negligible 
fraction of the coupling. For example, the standard recombination 
calculation yields a global ionized fraction Xi ^ 2 x_ 10"'' and 
Tk ~ 9[(1 + z)/20]^-*^ K for redshifts z ~ 10-100 dSeager et all 
1 19991) ; the exponent is slightly smaller than two (which would be 
expected for an adiabatically-cooling non-relativistic gas) because 
of a small amount of Compton heating off the CMB. At the higher 
redshifts, where Tk ~ 100 K, protons account for only ~ 0.04 per 
cent of the total coupling, but if the gas does remain cool to z ~ 20 
they provide ~ 2 per cent as much coupling as hydrogen-hydrogen 
collisions. In the unlikely event that the gas cools to T = 2.5 K 
at 2: = 10 without interference from luminous sources, we would 
have x^^/x™ ~ 0.6 - and protons would become even more im- 
portant than electrons. 

In reality, the first galaxies probably appear at 2; > 20. They 
flood the Universe with Lya photons (which affect the spin temper- 
ature through the Wouthuysen-Field mechanism) and X-ray pho- 
tons (which heat the gas). As a resu lt, the spin t emperature prob- 
ably increases well before z = 10 jSethill2005l ; lFurlanettdr2006l ; 
IPritchard & Furlanetto"2006^. so that protons never actually come 
to dominate the coupling. However, they still must be included for 
high-accuracy predictions of the 21 cm signal. 

In this paper, we have only examined spin-exchange via elec- 
tronically elastic proton-hydrogen collisions. At high temperatures, 
such collisions can instead excite higher atomic levels or even 
ionize the atom. Such processes could also affect the spin tem- 
perature because the hydrogen atom could enter a different spin 
state after the radiative cascade that follows. Although we do not 
know the relevant excitation cross sections, our experience with 
electron-hydrogen collisions suggests that in practice these colli- 
sions will not be important. For H-e", the most likely transition at 
relatively low temperatures is excitation to the 2p state, which is 
followed almost immediately by radiative de-excitation and emis- 
sion of a Lya photon. This Lya photon scatters ~ 10^ times be- 
fore redshifting out of resonance, so in terms of spin coupling it is 
much more important than the single collision that generated it. For 
a Maxwell-Boltzmann energy distribution, excitation of this pro- 
cess through electron-hydrogen collisions suddenly comes to dom- 
inate at Tk ~ 6300 K (Furlanetto & Furlanetto 2007). If proton- 
hydrogen collisions cause similar excitations, they too would be- 
come relatively unimportant once electronic excitations become 
energetically feasible. Furthermore, because the electron-hydrogen 
rate coefficient exceeds the proton-hydrogen rate coefficient by a 
factor of several in this temperature range, we expect proton colli- 
sions to be only a minor perturbation. At higher temperatures, the 
coUisionally-generated Lya photons completely dominate, and the 
effect of will be even smaller^ 

We thank P. Krstic for making his H-H+ cross section data 
available to us in an electronic form. This publication has been ap- 
proved for release as LA-UR-07-0513. Los Alamos National Lab- 
oratory, an affirmative action/equal opportunity employer, is oper- 
ated by Los Alamos National Security, LLC, for the National Nu- 
clear Security Administration of the U.S. under contract DE-AC52- 
06NA25396. 



Note that the decline in /k^^ atT> 3000 K is somewhat artificial. 
iFurla netto & Furlanetto (2007) excluded all electrons with E > 10.2 eV 
from the thermal average, because these electrons have a substantial proba- 
bility of exciting the atom to higher electronic states. 



In reality. X-ray heating of the IGM results in an excess of fast electrons 
even when the temperature is much smalle r. Thus the Lya channel can be 
important throughou t the "reheating" era Chen & Miralda-Escudd |2006| ; 
iChuzhov et ali20(33 ; |Pritchard & Furlanettdl2006lx 
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